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Abstract 

It  is  generally  accepted  that  the  ideal  operating  temperature  of  a  molten  carbonate  fuel  cell  (MCFC)  is  650  °C.  Nevertheless,  when  waste 
heat  utilization  in  the  form  of  an  expander  and  steam  production  cycle  is  introduced  in  the  system,  another  temperature  level  might  prove 
more  productive.  This  article  is  a  first  attempt  to  the  optimization  of  MCFC  operating  temperatures  of  a  MCFC  system  by  presenting 
a  case  study  in  which  the  efficiency  of  a  combined  heat  and  power  (CHP)  plant  is  analyzed.  The  fuel  cell  plant  under  investigation  is 
designed  around  a  250kW-class  MCFC  fuelled  by  natural  gas,  which  is  externally  reformed  by  a  heat  exchange  reformer  (HER).  The 
operating  temperature  of  the  MCFC  is  varied  over  a  temperature  range  between  600  and  700  °C  while  keeping  the  rest  of  the  system  the 
same  as  far  as  possible.  Changes  in  energetic  efficiency  are  given  and  the  causes  of  these  changes  are  further  analyzed.  Furthermore,  the 
exergetic  efficiencies  of  the  system  and  the  distribution  of  exergy  losses  in  the  system  are  given.  Flowsheet  calculations  show  that  there  is 
little  dependency  on  the  temperature  in  the  first  order.  Both  the  net  electrical  performance  and  the  overall  exergetic  performance  show  a 
maximum  at  approximately  675  °C,  with  an  electrical  efficiency  of  51.9%  (LHV),  and  an  exergy  efficiency  of  58.7%.  The  overall  thermal 
efficiency  of  this  CHP  plant  increases  from  87. 1%  at  600  °C  to  88.9%  at  700  °C.  Overall,  the  change  in  performance  is  small  in  this  typical 
range  of  MCFC  operating  temperature. 

©  2003  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  play  an  important  role  in  the  continuing  effort 
to  increase  the  efficiency  of  electricity  production  and  to  re¬ 
duce  atmospheric  pollution.  High  net  power  efficiency  can 
be  achieved  thanks  to  the  principle  of  direct  conversion  of 
chemical  energy  to  electrical  energy,  and  thereby  avoiding 
the  extra  steps  of  combustion,  heat  transfer,  expansion  and 
generation  as  in  a  conventional  plant.  Furthermore,  when 
high-temperature  cells  like  the  molten  carbonate  fuel  cell 
(MCFC)  are  used,  additional  increase  in  the  overall  effi¬ 
ciency  can  be  obtained  by  proficient  residual  heat  utilization. 
Then,  in  order  for  the  total  efficiency  to  be  as  high  as  pos¬ 
sible,  it  is  possible  that  ideal  fuel  cell  operating  conditions 
will  not  coincide  with  the  optimal  total  system  performance. 
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It  is  generally  accepted  that  the  ideal  operating  temper¬ 
ature  of  an  MCFC  is  650  °C.  For  current  state  of  the  art 
fuel  cell,  this  temperature  is  the  best  compromise  between 
performance  and  endurance  [1-3].  Optimization  studies  of 
MCFC  system  efficiencies  are  therefore  done  in  the  fields 
of  cell  and  stack  configuration  [4-6]  and  system  configu¬ 
rations  [7,8].  However,  when  waste  heat  utilization  in  the 
form  of  an  expander  and  steam  generation  is  introduced 
in  the  system,  another  temperature  level  might  prove  to  be 
more  productive.  Previously,  we  performed  a  theoretical 
study  to  the  efficiencies  of  fuel  cell  systems  using  simple 
assumptions  for  the  irreversible  losses  [9,10].  The  results 
were  obtained  by  assuming  temperature  independent  losses 
for  the  recovery  of  the  exergy  of  the  residual  heat.  In  prac¬ 
tice,  the  efficiency  of  waste  heat  recovery  depends  on  the 
temperature  of  the  system.  Furthermore,  the  interactions 
between  the  fuel  cell  stack,  auxiliary  equipment  and  waste 
heat  recovery  system  were  previously  neglected.  Here,  we 
present  more  detailed  flowsheet  calculations  of  a  typical 
external  reformed  MCFC  combined  heat  and  power  (CHP) 
plant  with  which  we  investigate  the  influence  of  the  fuel 
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Nomenclature 

Acell 

active  cell  area  (m2) 

Q 

fitting  constant  for  quasi-ohmic 

resistance  (£2  m2) 

d 

electrolyte  thickness  (mm) 

do 

normalization  constant  for  the 

electrolyte  thickness  (mm) 

Ex 

exergy  (kW) 

A  hi 

activation  enthalpy  (J/mol) 

feel  1 

current  load  density  of  the  unit  cell  (A/m2) 

m 

molar  fraction  (-) 

P 

power  delivered  (kW) 

P 

pressure  (bar) 

Pi 

partial  pressure  of  gas  of  species  i  (bar) 

A p 

pressure  loss  (bar) 

r 

quasi-ohmic  resistance  (£2  m2) 

R 

universal  gas  constant  (J/mol  K) 

T 

temperature  (°C) 

To 

standard  ambient  temperature  (°C) 

A  Thigh 

high  end  temperature  difference  (K) 

A7jow 

low  end  temperature  difference  (K) 

u 

local  cumulative  fuel  utilization  (-) 

Uf 

total  fuel  utilization  (-) 

Fcell 

cell  voltage  (V) 

Veq 

Nernst  voltage  (V) 

Greek  symbols 

0 

mass  flow  (kg/s) 

T] 

efficiency  (-) 

k 

air  to  fuel  ratio  (kg/kg) 

cell  operating  temperature  on  the  system  performance.  This 
system  has  been  modeled  and  implemented  in  the  program 
Cycle-Tempo  [11],  which  is  a  program  that  Delft  Univer¬ 
sity  has  developed  for  flowsheet  calculations.  The  effects 
of  varying  the  cell  temperature  upon  the  fuel  cell  (both  re¬ 
versible  and  kinetical)  and  various  system  aspects  have  been 
investigated,  and  based  on  these  results,  the  optimal  operat¬ 
ing  temperature  for  maximum  efficiency  will  be  presented. 


2.  System  configuration 

The  system  selected  for  this  study  is  similar  to  a 
system-design  considered  for  a  250  kW  natural  gas  MCFC 
system  as  jointly  defined  in  the  past  by  the  Delft  Univer¬ 
sity  of  Technology  and  ECN  (Netherlands  Energy  Research 
Foundation)  [12].  It  has  the  following  main  features: 

•  250  kW-class  CHP  system; 

•  natural  gas  as  primary  fuel; 

•  fuel  gas  is  externally  reformed; 

•  pressurized  system  operating  at  4  bar. 


The  initially  proposed  system  has  an  intricate  connection 
between  the  anode  cycle  and  the  cathode  cycle  (via  the  pre¬ 
heating  stages  and  mixing).  This  connection  is  removed  here 
mainly  as  regards  to  stability,  constructive  simplicity  and 
controllability.  Separating  the  two  flows  should  bring  about 
a  better  system.  Fig.  1  shows  the  flowsheet  of  the  mod¬ 
ified  system  that  is  taken  directly1  from  the  flowsheeting 
program.  The  system  layout  shows  that  apart  from  the  fuel 
cell  we  can  distinguish  five  subsystems.  Next,  the  fuel  cell 
and  the  subsystems  will  be  introduced  by  explaining  their 
functions.  Furthermore,  the  input  data  used  for  the  different 
components  will  be  presented.  These  inputs  consist  of  input 
data  that  characterize  state  of  the  art  equipments. 

2.7.  Fuel  cell 

A  unique  feature  of  the  fuel  cell  model  in  the  flowsheeting 
program  is  its  capability  of  calculating  design  and  off  de¬ 
sign  or  part  load  performances,  next  to  the  energy  and  mass 
transfer  calculations.  The  cell  performance  is  calculated  by 
numerically  solving  an  integral  expression  for  cell  voltage 
Vceii  as  function  of  the  operating  parameters  fuel  utilization 
Uf  and  current  density  /ceii 

1  fu f  r  p  , 

Vceii  =  —  /  veq(u)du-- — /  i2(x)  dx  (1) 

Uf  JO  Aell  Jo 

where  Veq  is  the  Nernst  voltage  expressed  as  function  of 
the  cumulative  local  fuel  utilization  u ,  i(x)  the  local  current 
density  and  r  the  quasi-ohmic  resistance  that  accounts  for 
all  irreversible  losses.  Therefore,  the  first  integral  represents 
the  reversible  cell  voltage  as  function  gas  composition,  op¬ 
erating  temperature  and  pressure  and  total  fuel  utilization 
Uf.  The  second  integral  represents  the  irreversible  losses 
due  to  ohmic  losses  and  electrochemical  kinetics  as  func¬ 
tion  of  mean  current  density  /ceii .  This  is  done  by  assuming 
a  local  ohmic  relation  for  the  electrode  kinetics  (see  [13]). 
A  detailed  description  of  the  complete  fuel  cell  model  is 
given  in  [14]  together  with  a  verification  of  model  by  com¬ 
paring  the  calculated  results  with  the  experimental  results 
obtained  from  a  110  cm2  benchmarking  class  MCFC  single 
cell.  Using  the  measured  macroscopic  cell  resistance  as  the 
quasi-ohmic  resistance,  this  verification  showed  an  average 
relative  discrepancy  of  0.5%  over  a  wide  range  of  operating 
condition,  and  a  maximum  discrepancy  of  3%  at  full  load. 
Hence,  it  is  shown  the  fuel  cell  model  is  correct  and  its 
accuracy  sufficient  to  be  used  for  flowsheeting  purposes. 
This  fuel  cell  model  enables  us  to  use  the  relations  for  the 


1  Components  that  are  needed  solely  for  starting  values  of  mathematical 
iterations  are  removed  here. 

2  A  small  difference  in  value  exists  between  the  quasi-ohmic  resistance 
and  the  macroscopic  cell  resistance.  The  difference  is  described  in  detail 
in  [13].  Using  the  fitted  quasi-ohmic  resistance,  the  discrepancy  between 
the  calculated  and  measured  cell  voltage  can  be  improved  to  1.7%  at 
full  load  (see  [14]).  The  elaborate  fitting  procedure  for  the  quasi-ohmic 
resistance  is  here  omitted. 
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temperature  dependency  of  the  macroscopic  cell  resistance 
measured  by  CRIEPP  [6,15].  They  obtained  the  following 
empiric  relations  for  the  cell  resistances  r,  which  can  be 
distinguished  in  an  anode  contribution  ra,  an  electrolyte 
resistance  rlv  and  a  cathode  contribution  rc : 

r  =  ra  +  rir  +  rc  (2) 


with: 


r  a  =  CaT  exp 


1 

y/Pz,  H2  ’ 


r\ 


ir 


rc  =  C ci  exp 


AAd 

RT 


T  VPco2 

0.75 

Pc,  02 


CC2  exp(A/zC2/R7)^ 

^c,C02  +  Cc3  exp(A/zc3/^7)^c,H20 


The  symbols  Q,  A/z/  and  Jb  are  respectively  fitting  para¬ 
meters,  activation  enthalpy  and  normalization  parameter 
for  the  electrolyte  matrix  thickness  with  the  thickness  d  = 


3  Central  Research  Institute  of  Electric  Power  Industry,  Japan. 


0.916  mm.  The  fitting  parameters  and  activation  enthalpy 
are  obtained  by  fitting  these  empirical  relations  to  exper¬ 
imental  measurements  performed  at  a  temperature  range 
between  600  and  700  °C,  and  at  a  pressure  range  between 
1  and  5  bar.  The  resulting  values  are  given  in  Table  1 . 

The  empirical  relations  are  used  to  calculate  the  cell  resis¬ 
tance  as  function  of  cell  temperature,  pressure  and  average 


Table  1 


Values  for  the  fitting  parameters  for  the  cell  resistance  taken  from  [6,14] 


Parameter 

Value 

Ca 

9.50  x  10-7  £2  cm2 

Cel 

6.91  x  10-15  £2  cm2 

Cc2 

3.75  x  10~9  Q  cm2 

CC3 

1.07  x  10~6  Q  cm2 

Cir 

9.48  x  10“ 3  £2  cm2 

Aha 

27.9kJ/mol 

A/zci 

179.2  kJ/mol 

Ahc2 

67.2kJ/mol 

Ahc3 

95.2  kJ/mol 

A/lir 

23.8  kJ/mol 

dfo 

0.916mm 
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Table  2 


Input  parameters  of  the  system 


Fuel  cell 

^tcell 

250  m2 

Cell 

1500  A/m2 

Uf 

70% 

Tide— ac 

96% 

P 

4  bar 

A/?  anode 

0.05  bar 

A/?Cathode 

0. 1  bar 

Tout  Tn 

100  °C 

Anode  gas  recirculation  and  moisture  separation 

Tout  anode  recirculation  gas 

460  °C 

Thot  water 

80°C 

Tutilization 

20  °C 

A  Flow  evaporator 

20  °C 

A p  heat  exchangers 

0.05-0.10  bar 

A p  moist  separator 

Primary  side  0.15  bar; 
secondary  side  0.1  bar 

iji  for  pump 

0.70 

rji  for  blower 

0.72 

Heat  exchange  reformer  and  fuel  preheat 

k  (air  factor)  combustor 

1.1 

Steam  to  fuel  ratio  reformer 

2.59 

Treact  reformer 

800  °C 

/? react  reformer 

4  bar 

A p  reformer 

Primary  side  0.5  bar; 
secondary  side  0.25  bar 

Tout  reformer 

According  to  cell 
inlet  temperature 

A p  heat  exchangers 

0.05  bar  for  both  primary 
and  secondary  sides 

Tfeed  reformer 

480  °C 

Fair  combustor 

440  °C 

r]i  compressor 

0.72 

Cathode  gas  recirculation 

A p  heat  exchanger 

0.05  bar  for  both  primary 
and  secondary  sides 

r]i  compressor 

0.72 

Expander  and  waste  heat  boiler 

A  Flow  evaporator 

20  °C 

A p  heat  exchanger 

Primary  side  0.1  bar; 
secondary  side  0.01  bar 

A p  utilization 

0.2  bar 

r]i  expander 

0.75 

rj  generator 

0.9 

r]i  pumps 

0.75 

gas  composition  (by  means  of  average  partial  pressures  pi 
and  mol  fractions  mf).  The  cell  resistance  r  determines  the 
irreversible  losses  and  therefore  the  performance  of  the  cell. 
Changes  in  reversible  heat  production  and  Nernst  loss  due 
to  temperature  and  gas  compositions  changes  are  accounted 
by  the  first  integral  in  Eq.  (1)  by  means  of  changes  in  the 
local  Nernst  voltages  [13].  Note  that  the  symbols  p  and  T  in 
Eq.  (3)  are  used  for  respectively  partial  pressure  (bar)  and 
absolute  temperature  (K). 

In  this  study,  both  the  cell  area  Acen  and  the  current  den¬ 
sity  /cell  are  kept  constant  and  they  are  given  in  Table  2. 
The  power  that  is  delivered  by  the  fuel  cell  is  therefore 


a  function  of  cell  resistance,  fuel  utilization  and  gas  inlet 
compositions.  Losses  due  to  the  dc  to  ac  conversion  are  in¬ 
troduced  by  the  inverter  efficiency  77dc-ac-  Fuel  utilization  is 
fixed  at  70%,  which  is  5%  lower  than  described  in  [12].  The 
reason  for  the  low  fuel  utilization  lays  solely  on  the  sepa¬ 
ration  of  anode  and  cathode  cycles.  Due  to  this  separation, 
extra  enthalpy  is  required  in  order  to  heat  up  the  fuel  and  to 
provide  enough  heat  for  the  heat  exchange  reformer  (HER) 
without  subtracting  this  from  the  cathode  cycle.  By  reducing 
the  fuel  utilization,  the  fuel  input  has  increased  and  more 
heat  is  available  from  the  anode  off  gas.  Friction  losses  are 
introduced  by  imposing  pressure  drops  of  0.05  and  0. 1  bar 
for  the  anode  and  cathode  respectively.  Co-flow  configura¬ 
tion  is  assumed  and  the  temperature  difference  between  the 
inlets  and  the  outlets  of  both  the  anode  and  cathode  are  set 
at  100  K. 

The  system  is  implemented  in  such  a  way  that  the  fuel 
cell  is  the  dominant  apparatus  and  that  both  fuel  and  air  con¬ 
sumptions  are  mainly  determined  by  it.  Fuel  consumption 
is  determined  by  /ce n,  Aceu  and  Uf.  The  cathode  mass  flow, 
and  the  related  air  consumption,  is  determined  by  the  heat 
that  must  be  discharged  from  the  fuel  cell. 

2.2.  Anode  gas  recirculation  and  moisture  separation 

Fuel  that  is  not  converted  by  the  fuel  cell  is  combusted 
in  the  reformer  (HER).  However,  the  anode  off  gas  con¬ 
tains  large  amount  of  moisture  that  will  adversely  influence 
the  performance  of  the  HER.  The  anode  off  gas  is  therefore 
cooled  in  several  stages  to  separate  most  of  the  moisture.  The 
transferred  heat  is  used  for  heating  up  and  evaporating  wa¬ 
ter  that  is  needed  for  the  reforming  reaction.  Heat  released 
in  the  moisture  separator  is  utilized  by  external  consumers 
(e.g.  a  district  heating  system)  represented  here  by  a  heat 
sink.  The  moisture  separator  produces  hot  water  at  80  °C. 
After  utilization,  this  water  returns  at  a  temperature  of  60  °C 
and  it  is  recirculated  back  to  the  moisture  separator.  The  an¬ 
ode  off  gas  is  circulated  by  a  blower,  reheated  and  sent  to 
the  HER.  This  dried  anode  recycle  gas  leaves  the  anode  gas 
recirculation  and  moisture  separation  subsystem  at  a  fixed 
temperature  of  460  °C  in  order  to  keep  the  inlet  tempera¬ 
ture  of  the  heat  exchange  reformer  constant  throughout  this 
study. 

The  efficiency  of  the  heat  transfer  processes  depends 
strongly  on  the  choice  of  flow  configuration  and  the  final 
temperature  differences  between  primary  and  secondary 
flows.  Here,  all  heat  exchangers  are  operated  in  counter  flow 
mode.  Only  the  low  end  temperature  difference  (A7iow)  of 
the  steam  evaporator  (apparatus  2)  is  set  at  20  K.  The  others 
are  calculated  using  the  fixed  inlet  temperature  of  the  HER 


4  See  [15]  for  more  details  of  fuel  cell  modeling. 

5  The  other  consumer  of  air  is  the  combustion  chamber  of  the  heat 
exchange  reformer.  This  amount  of  air  is  relatively  small. 
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and  the  boiling  temperature  of  water  at  the  exchanger’s  out¬ 
let. 

2.3.  Heat  exchange  reformer  and  fuel  preheat 

The  heat  exchange  reformer  is  modeled  here  by  a  combus¬ 
tion  chamber  and  a  steam-reforming  reactor.  The  combus¬ 
tion  chamber  is  fuelled  by  the  dried  anode  off  gas  and  the  air 
factor  X  of  combustion  is  set  at  1 . 1  with  which  flue  gas  at  a 
temperature  of  over  1250  °C  is  obtained.  The  heat  that  can 
be  derived  from  the  flue  gas  is  sufficient  for  the  reforming 
reaction  (apparatus  9,  Fig.  1),  superheating  steam  (apparatus 
18)  and  preheating  fuel  (apparatus  19).  The  remaining  heat 
is  used  for  heating  air  (apparatus  17)  that  is  supplied  to  the 
combustion  chamber  (apparatus  10).  The  temperature  of  the 
air  leaving  the  heat  exchanger  is  set  at  440  °C  and  the  tem¬ 
perature  of  both  natural  gas  and  steam  that  enter  the  HER 
are  set  at  480  °C.  The  natural  gas  heater  (apparatus  19)  and 
the  steam  heater  (apparatus  18)  are  placed  here  in  parallel 
since  both  heat  exchangers  are  in  practice  combined  in  a 
single  unit. 

After  passing  the  air  preheater  (apparatus  17),  the  CO2 
rich  flue  gas  is  mixed  with  the  recycled  cathode  gas  and 
preheated  fresh  air  to  provide  the  MCFC  cathode  with  O2 
and  CO2. 

The  reforming  reaction  is  modeled  by  assuming  chemical 
equilibrium  at  800  °C  and  4  bar.  The  ratio  of  steam  to  fuel 
is  set  here  to  2.59  kg/kg.  Friction  losses  are  introduced  in 
the  reformer  by  imposing  pressure  drops  of  0.5  and  0.25  bar 
for  the  primary  process  flow  (the  product  gas  flow)  and  the 
secondary  heat  exchange  flow  (the  flue  gas  flow)  respec¬ 
tively.  Other  friction  losses  are  introduced  by  pressure  drops 
of  0.05  bar  for  both  primary  and  secondary  sides  of  the  fuel 
preheating  line.  The  isentropic  efficiency  of  the  air  compres¬ 
sor  is  assumed  to  be  0.72. 

2.4.  Cathode  gas  recirculation 

The  cathode  gas  not  only  provides  O2  and  CO2  for  the 
electrochemical  reaction,  it  also  serves  as  the  main  coolant 
for  the  fuel  cell,  and  therefore  the  mass  flow  of  the  cath¬ 
ode  gas  has  to  meet  the  cooling  requirements.  This  mass 
flow  of  air  necessary  for  cooling  is  far  greater  than  required 
for  the  cathode  reaction.  Part  of  this  air  is  therefore  recircu¬ 
lated  and  the  amount  of  recirculation  is  set  accordingly  to 
assure  the  fixed  cathode  inlet  temperature  after  mixing  this 
recirculation  flow  with  fresh  air  and  flue  gas  from  the  HER. 
Before  mixing,  this  recycle  flow  is  partly  cooled  by  preheat¬ 
ing  the  pressurized  fresh  air.  The  flue  gas  from  the  HER 
is  the  main  source  of  CO2  required  for  the  cathode  reac¬ 
tion.  In  all  considered  situations,  the  concentration  of  CO2  at 
the  cathode  inlet  is  above  the  commonly  assumed  minimum 
of  8  mol%. 


6  This  means  10%  more  air  than  needed  for  stochiometric  combustion. 


2.5.  Expander  and  waste  heat  boiler 

The  hot  and  pressurized  gas  that  leaves  the  cathode 
recycle  loop  produces  electricity  through  an  expander  and 
the  attached  generator.  Losses  are  introduced  by  defining 
isentropic  efficiency  for  the  expander  (here  75%)  and  con¬ 
version  efficiency  for  the  electrical  generator  (here  95%). 
After  expansion,  the  temperature  of  the  flue  gas  is  suffi¬ 
ciently  high  to  produce  saturated  steam.  The  pressure  and 
temperature  of  this  steam  is  set  at  lObar  and  about  180  °C. 
This  steam  can  be  applied  for  industrial  heating  purposes 
and  the  utilization  of  this  heat  is  represented  here  by  a  heat 
sink.  The  returning  condensate  from  this  sink  is  used  to 
feed  the  boiler.  Finally,  the  residual  flue  gas  is  discharged 
to  the  environment  via  a  flue  gas  stack. 

3.  Input  data  and  calculations 

The  system  performance  depends  strongly  on  the  in¬ 
put  data.  Especially  isentropic  efficiencies  of  rotating 
equipment,  pressure  drops  and  pinch  points  of  heat  ex¬ 
changers  determine  the  irreversible  losses  of  the  system 
and  therefore  also  the  calculated  overall  efficiencies.  For 
this  study,  we  have  used  a  combination  of  input  data 
that  characterize  state  of  the  art  equipment.  Table  2  gives 
the  main  input  parameters  for  the  components  in  the 
system. 

The  energy  input  to  the  system  is  determined  by  the  size 
of  the  fuel  cell,  anode  gas  composition  and  fuel  utiliza¬ 
tion.  Since  these  parameters  are  fixed  for  all  calculations, 
the  energy  input  of  the  system  is  constant  and  in  all  cases 
the  energy  input  is  557.57  kW,  based  on  the  lower  heat¬ 
ing  value  (LHV),  and  exergy  input  is  580.82  kW  based  on 
7q  =  25  °C.  The  source  for  the  energy  and  exergy  input 
is  the  natural  gas,  which  is  of  Dutch  “Slochteren”  quality 
with  as  main  components  about  81  mol%  CH4  and  14mol% 
N2  and  with  a  LHV  of  708.22  kJ/mol.  Other  mass  input  of 
this  plant  is  air  which  is  defined  according  to  the  ISO  stan¬ 
dard.  The  exact  compositions  that  are  used  for  both  natural 
gas  and  air  can  be  found  in  the  handbook  of  the  program 
[11]. 

The  fuel  cell  system  is  analyzed  at  five  different  cell  oper¬ 
ating  temperatures,  i.e.  600,  625,  650,  675  and  700  °C.  Using 
650  °C  as  a  reference,  the  operating  temperature  of  the  fuel 
cell  is  adjusted  by  changing  the  cell  temperature  and  the 
quasi-ohmic  resistance  of  the  cell  together  with  the  change  in 
the  cathode  gas  recirculation  percentage.  Any  change  in  the 
fuel  cell  temperature  influences  the  quasi-ohmic  resistance, 
which  on  its  turn  influences  both  cathode  gas  flow  and  the 
amount  of  the  cathode  gas  recirculation.  Consequently,  the 
recirculation  influences  the  cathode  gas  composition,  which 
determines  the  quasi-ohmic  resistance.  All  these  parameters 
are  therefore  closely  related  and  several  manually  controlled 
iterations  are  needed  in  order  to  find  the  solution  for  each 
temperature. 
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4.  Results  and  discussion 

Starting  with  the  overall  results,  Fig.  2  shows  the  overall 
system  thermal  efficiency  (77totalLHv)  and  output  distribu¬ 
tion  based  on  the  energy  input.  Note  that  the  surfaces  in  this 
scale  do  not  represent  the  ratio  of  power  over  heat  correctly 
since  we  adapted  the  scale  to  emphasize  different  results. 
The  numerical  values  are  also  summarized  in  Table  3.  Here, 
we  should  note  that  all  numbers  in  the  tables  are  given  in 
at  least  two  digits  behind  the  decimal  point.  This  suggests  a 
high  level  of  precision  in  our  computer  simulations.  On  the 
other  hand,  we  have  used  several  estimated  input  values  for 
the  performance  of  heat  exchangers  and  rotating  equipment 
and  consequently  the  absolute  precision  in  the  calculated 

Table  3 


Summary  of  energy  output 


T;ell 

600  °C 

625  °C 

650  °C 

675  °C 

700  °C 

Pfc  (kW) 

297.54 

304.20 

306.19 

305.95 

303.86 

^expander  (kW) 

73.86 

65.60 

60.05 

58.64 

58.20 

P aux  (kW) 

-93.98 

-84.08 

-77.58 

-75.28 

-74.22 

Pm t  (kW) 

277.42 

285.72 

288.66 

289.32 

287.84 

7netLHV  (%) 

49.76 

51.24 

51.77 

51.89 

51.62 

P steam  (kW) 

100.24 

94.67 

91.59 

93.69 

96.93 

P water  (kW) 

107.86 

108.68 

109.49 

110.28 

111.06 

Aieat  (kW) 

208.09 

203.35 

201.08 

203.97 

207.99 

7?heat LHV  (%) 

37.32 

36.47 

36.06 

36.58 

37.30 

P total  (kW) 

485.52 

489.07 

489.74 

493.29 

495.82 

7 total  LHV  (%) 

87.08 

87.72 

87.84 

88.47 

88.93 

efficiencies  can  therefore  not  be  guaranteed  in  practice.  Nev¬ 
ertheless,  these  numbers  are  not  round  off  in  further  extend 
since  otherwise  the  difference  in  the  calculated  results  will 
not  become  apparent.  Since  all  calculations  are  based  on  the 
same  system  using  the  same  system  inputs,  the  qualitative 
result  is  therefore  not  affected  by  the  estimated  inputs. 

First,  we  see  that  the  electrical  efficiency  (77netLHv)  of  the 
system  increases  with  operating  temperature  and  reaches 
a  maximum  at  675  °C.  The  difference  between  the  maxi¬ 
mum  and  minimum  efficiency  is  about  2.1%  point.  Second, 
the  overall  efficiency  (i.e.  heat  and  power)  increases  with 
operating  temperature  as  well  but  we  did  not  find  a  peak 
maximum  in  the  temperature  range  we  investigated.  In  this 
temperature  range,  the  difference  between  the  maximum 
and  minimum  is  here  about  1.9%  point.  Finally,  Table  3 
shows  that  the  auxiliary  power  consumption  Pmx  decreases 
with  operating  temperature  from  about  25%  of  the  gross 
power  production  at  600  °C  to  about  20%  at  700  °C.  The 
auxiliary  power  is  mainly  used  by  the  compressor  for  com¬ 
pressing  fresh  air  for  the  cathode,  and  the  change  in  this 
causes  the  differences  in  auxiliary  power  consumption.  The 
rest  of  the  auxiliary  power  is  mainly  used  by  the  compressor 
for  air  to  feed  the  HER  and  by  the  two  recycle  blowers. 

Next,  the  causes  of  the  changes  in  fuel  cell  power  output 
are  discussed  in  more  detail.  Fig.  3  shows  the  cell  resis¬ 
tance  and  the  net  power  delivered  by  the  fuel  cell  stack.  As 
expected,  the  irreversible  losses  given  by  the  cell  resistance 
decreases  with  increasing  cell  temperature.  This  results  to 
an  increase  in  stack  performance  and  thus  an  increase  in 
delivered  power.  On  the  other  hand,  the  reversible  open 
cell  voltage  (OCV)  given  by  the  Nernst  equation  decreases 
linearly  with  increasing  temperature,  as  described  in  details 
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in  our  previous  theoretical  study  [9].  The  opposite  temper¬ 
ature  behavior  of  irreversible  losses  and  OCV  results  here 
to  a  maximum  in  power  output  at  rceii  =  650  °C.  In  our 
previous  theoretical  study  [9],  we  did  not  find  a  maximum 
in  cell  performance  in  this  typical  range  of  operating  tem¬ 
perature.  There,  we  used  a  constant  gas  composition  while 
in  this  present  study  the  cathode  gas  composition  is  mainly 
determined  by  the  cooling  requirement  of  the  stack  and 
by  the  heat  requirement  of  the  cathode  gas  recirculation 
subsystem.  Therefore,  the  exact  cathode  gas  composition  is 
here  a  function  of  operating  temperature  and  this  gas  com¬ 
position  influences  both  the  cell  resistance  as  well  as  the 
OCV.  The  difference  in  the  optimum  in  stack  performance 
of  our  present  results  and  our  previous  theoretical  result  can 
be  ascribed  to  the  differences  in  cathode  gas  composition. 

The  flowsheet  calculations  show  that  the  overall  system 
has  a  different  optimum  operating  temperature  than  the  fuel 
cell  stack.  This  is  caused  by  the  net  expander  output  ^expander 
and  auxiliary  power  consumption  Pmx.  We  have  analyzed 
this  by  examining  the  cell  resistance  and  cathode  recircula¬ 
tion  data,  both  given  in  Table  4.  First,  we  note  that  the  dif¬ 
ference  in  net  performance  (77netLHv),  in  particular  between 
650  and  700  °C  is  very  small  (a  difference  of  only  0.27%). 


This  result  is  in  accordance  with  our  observation  for  the 
whole  system  as  considered  in  that  same  theoretical  study 
[9].  First,  at  low  temperature,  the  net  performance  is  ad¬ 
versely  affected  by  the  high  irreversible  losses  while  at  high 
operating  temperature  the  system  performance  remains  con¬ 
stant.  Second,  high  electrical  output  by  the  fuel  cell  should 
result  to  low  heat  release  and  consequently  little  cooling  is 
required.  From  600  to  650  °C,  Table  4  shows  a  decreasing 
^cathode*  which  is  a  direct  result  of  the  increasing  Ppc-  It  is 
however  interesting  to  note  that  although  the  electrical  output 
Fpc  is  highest  at  650  °C,  the  cathode  mass  flow  ^cathode  and 
cooling  requirement  is  lowest  at  675  °C.  This  seems  to  con¬ 
tradict  what  is  expected  from  theory,  since  the  highest  power 
output  should  result  in  the  lowest  cooling  requirement.  The 
cause  of  this  contradiction  can  be  ascribed  to  the  difference 
in  gas  composition  in  the  anode  outlet  due  to  difference  in 
equilibrium  in  the  hydrogen- shift  reaction.  At  650  °C  the 
average  enthalpy  of  the  anode  outlet  is  33.31  kJ/mol  (LHV) 
while  at  675  °C  it  is  33.42  kJ/mol  (LHV).  Since  the  anode 
mass  flow  is  constant  for  all  calculations,  the  enthalpy  re¬ 
lease  by  anode  outlet  is  slightly  higher  675  °C.  The  higher 
enthalpy  release  reduces  the  cooling  requirement  of  the  fuel 
cell  and  hence  lower  cathode  mass  flow.  The  difference  in 


Table  4 


Summary  of  quasi-ohmic  resistance  and  cathode  recirculation 


Tell 

600  °C 

625  °C 

650  °C 

675  °C 

700  °C 

r  (£2  cm2) 

0.9619 

0.7501 

0.6072 

0.4976 

0.4187 

^recirculation 

72.83 

75.02 

77.09 

78.00 

78.84 

^cathode  (kg/s) 

1.383 

1.304 

1.267 

1.256 

1.262 

expander  (kg/ s) 

0.344 

0.297 

0.264 

0.251 

0.243 

Table  5 

Summary  of  exergy  output  (with  Tq  =  25  °C) 


Tell 

600  °C 

625  °C 

650  °C 

675  °C 

700  °C 

EXnet  (kW) 

277.42 

285.72 

288.66 

289.32 

287.84 

EXsteam  (kW) 

36.48 

34.46 

33.34 

34.10 

35.28 

Exwater  (kW) 

17.27 

17.41 

17.54 

17.66 

17.79 

7Ex  (%) 

57.02 

58.12 

58.46 

58.72 

58.69 
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cathode  mass  flow  turns  the  optimum  temperature  for  the 
system  to  675  °C  from  the  optimum  temperature  of  650  °C 
for  the  stack.  Finally,  this  study  also  shows  that  operating 
at  elevated  temperature  requires  increase  in  cathode  recy¬ 
cling  due  to  the  higher  inlet  temperature  of  the  fuel  cell. 
This  is  shown  in  Table  4  where  the  recirculation  percentage 


Of  the  Cathode  gas  %recirculation,  Cathode  mass  flow  ^cathode 
and  expander  mass  flow  Expander  are  given.  The  increase 
in  %  recirculation  reduces  expander  power  output  and  the  over¬ 
all  auxiliary  power  consumption.  The  latter  is  due  to  the 
lower  air  input  and  less  work  required  for  the  compression 
of  air. 


Fig.  5.  Exergy  loss  of  subsystems  as  function  of  temperature. 
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Table  3  gives  the  overall  results  and  it  lists  the  amount  of 
useful  heat  produced  by  the  system  that  is  consumed  by  ex¬ 
ternal  users.  The  production  of  hot  water  at  80  °C  increases 
monotone  with  operating  temperature  while  the  steam  pro¬ 
duction  at  180  °C  shows  a  minimum  at  650  °C.  The  change 
of  the  latter  dominates  the  overall  heat  production  resulting 
in  a  minimum  total  heat  production  (Pheat)  °f  201.08  kW  at 
650  °C.  The  overall  efficiency  ijtotai  (i.e.  combined  electric¬ 
ity  and  heat)  increases  monotone  with  operating  tempera¬ 
ture.  We  should  note  that  the  main  purpose  of  the  fuel  cell 
plant  is  the  production  of  electricity  while  the  production  of 
heat  is  of  minor  importance.  This  is  more  apparent  when  we 
evaluate  this  system  based  on  exergy.  Fig.  4  shows  the  exer- 
getic  efficiency  and  exergy  output  distribution  of  the  system 
(note  the  scale  of  this  figure).  The  numerical  values  are  sum¬ 
marized  in  Table  5.  It  is  clear  that  the  exergy  represented  by 
the  produced  heat  is  relatively  small.  Based  on  exergy,  the 
system  efficiency  r/EX  shows  a  maximum  at  675  °C.  The  dif¬ 
ference  in  overall  exergetic  efficiency  between  the  highest 
and  lowest  value  is  here  1.7%.  Again,  the  change  in  exergy 
efficiency  in  the  temperature  range  between  650  and  700  °C 
is  small. 

Furthermore,  we  have  analyzed  the  exergy  loss  of  the  sub¬ 
systems  and  their  temperature  dependency  and  this  is  shown 
by  Fig.  5.  It  shows  that  the  exergy  loss  of  the  subsystems 
HER,  fuel  preheat  and  expander,  waste  heat  boiler  are  lit¬ 
tle  temperature  dependant.  The  exergy  loss  of  the  fuel  cell 
and  cathode  recycling  decreases  while  the  anode  recycling 
increases  with  operating  temperature.  The  overall  result  is  a 
minimum  of  exergy  loss  at  675  °C  and  thus  an  optimum  in 
overall  exergy  efficiency. 


Finally,  Fig.  6  shows  the  exergy  loss  distribution  at  the 
typical  operating  temperature  of  Tce\\  =  650  °C.  It  shows 
that  the  HER  has  the  highest  contribution  to  the  total  exergy 
loss.  The  combustion  process  of  the  HER  is  the  main  cause 
of  the  exergy  loss  of  this  system  (about  65%)  while  the 
contribution  from  heat  transfer  to  the  reforming  reaction 
is  relatively  small  (the  remaining  35%).  Improvement  in 
system  efficiency  is  therefore  expected  when  changing  this 
external  reforming  configuration  to  internal  reforming.  This 
off  course  will  involve  different  MCFC  stack  technology  and 
eventually  will  affect  the  system  layout. 

5.  Conclusion 

The  influence  of  the  operating  temperature  of  the  fuel  cell 
on  the  overall  system  efficiency  is  small  in  the  operating 
range  between  650  and  700  °C.  This  result  is  in  accordance 
with  the  results  of  our  previous  study  [9].  This  fuel  cell 
system  performs  best  at  675  °C  with  a  net  electrical  effi¬ 
ciency  of  51.89%  point  (based  on  LHV).  This  is  the  main 
conclusion  since  the  production  of  electricity  is  the  objec¬ 
tive  while  the  production  of  heat  is  of  minor  importance. 
Production  of  heat  plays  a  role  when  this  system  is  inte¬ 
grated  to  industrial  processes  together  with  district  heating 
system.  Exergy  conservation  is  then  an  additional  require¬ 
ment  for  a  sustainable  society.  The  exergetic  efficiency  of 
this  system  is  at  maximum  as  well  at  675  °C  with  a  value 
of  58.72%.  The  overall  energetic  CHP  efficiency  based  on 
LHV  increases  with  operating  temperature  and  the  highest 
value  is  achieved  at  the  highest  temperature  considered. 
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The  exergy  loss  contributed  by  the  heat  exchange  reformer 
is  the  highest  of  all  subsystems  and  this  should  be  tackled 
case  of  further  system  optimization. 

The  conclusions  regarding  the  efficiencies  as  function  of 
operating  temperature  as  presented  here  are  restricted  to  this 
specific  plant  design.  Nevertheless,  the  present  study  has 
shown  the  complexity  of  a  fuel  cell  system.  The  refinement 
by  detail  flowsheet  analysis  as  presented  here  has  revealed 
interactions  between  subsystems  that  cannot  be  seen  oth¬ 
erwise.  Examples  are  the  different  optimum  temperature 
for  the  fuel  cell  and  overall  system  and  the  mismatch  be¬ 
tween  the  optimum  stack  temperature  and  the  minimum 
stack-cooling  requirement.  Furthermore,  we  have  shown 
the  complex  interactions  between  the  different  processes 
in  a  system.  We  have  seen  here  that  changes  in  cell  tem¬ 
perature  involve  the  following  changes:  (1)  reversible  heat 
production  and  irreversible  losses  of  the  fuel  cell;  (2)  cool¬ 
ing  requirement  of  the  cell  and  consequently  the  auxiliary 
power  consumption;  and  (3)  changes  in  recirculation  mass 
flows  due  to  the  changes  in  fuel  cell  inlet  temperature.  Their 
relation  can  only  be  made  visible  by  flowsheet  calculations. 
This  study  has  therefore  shown  the  importance  of  flowsheet 
calculation  during  the  evaluation  of  the  complete  fuel  cell 
plant  when  changing  process  parameters.  Simple  theoretical 
calculations  can  show  first  order  trends  but  detailed  flow¬ 
sheet  calculations  are  required  due  to  the  complex  behavior 
and  intricate  interactions  in  a  fuel  cell  plant. 
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